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ABSTRACT 

The  effects  of  surface  wind  speed  and  direction  gradients  on  midlatitude  surface  vorticity  and  divergence 
fields  associated  with  mesoscale  sea  surface  temperature  (SST)  variability  having  spatial  scales  of  100- 
1000  km  are  investigated  using  vector  wind  observations  from  the  SeaWinds  scatterometer  on  the  Quick 
Scatterometer  (QuikSCAT)  satellite  and  SST  from  the  Advanced  Microwave  Scanning  Radiometer  for  Earth 
Observing  System  (AMSR-E)  Aqua  satellite.  The  wind-SST  coupling  is  analyzed  over  the  period  June  2002- 
August  2008,  corresponding  to  the  first  6+  years  of  the  AMSR-E  mission.  Previous  studies  have  shown  that 
strong  wind  speed  gradients  develop  in  response  to  persistent  mesoscale  SST  features  associated  with  the 
Kuroshio  Extension,  Gulf  Stream,  South  Atlantic,  and  Agulhas  Return  Current  regions.  Midlatitude  SST 
fronts  also  significantly  modify  surface  wind  direction;  the  surface  wind  speed  and  direction  responses  to 
typical  SST  differences  of  about  2°-4°C  are,  on  average,  about  1-2  m  s~^  and  4°-8°,  respectively,  over  all  four 
regions.  Wind  speed  perturbations  are  positively  correlated  and  very  nearly  collocated  spatially  with  the  SST 
perturbations.  Wind  direction  perturbations,  however,  are  displaced  meridionally  from  the  SST  perturba¬ 
tions,  with  cyclonic  flow  poleward  of  warm  SST  and  anticyclonic  flow  poleward  of  cool  SST. 

Previous  observational  analyses  have  shown  that  small-scale  perturbations  in  the  surface  vorticity  and 
divergence  fields  are  related  linearly  to  the  crosswind  and  downwind  components  of  the  SST  gradient,  re¬ 
spectively.  When  the  vorticity  and  divergence  fields  are  analyzed  in  curvilinear  natural  coordinates,  the  wind 
speed  contributions  to  the  SST-induced  vorticity  and  divergence  depend  equally  on  the  crosswind  and 
downwind  SST  gradients,  respectively.  SST-induced  wind  direction  gradients  also  significantly  modify  the 
vorticity  and  divergence  fields,  weakening  the  vorticity  response  to  crosswind  SST  gradients  while  enhancing 
the  divergence  response  to  downwind  SST  gradients. 


1.  Introduction 

On  spatial  scales  spanning  ocean  basins,  sea  surface 
temperature  (SST)  perturbations  are  found  to  be  nega¬ 
tively  correlated  with  surface  wind  speed  perturbations 
(e.g.,  Mantua  et  al.  1997;  Okumura  et  al.  2001;  Xie  2004). 
Over  these  broad  spatial  scales,  large-scale  atmospheric 
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circulation  patterns  change  surface  ocean  temperatures 
through  modulation  of  surface  heat  fluxes  and  upper 
ocean  mixing  (e.g.,  Cayan  1992).  On  smaller  spatial 
scales  between  100  and  1000  km,  however,  surface  wind 
speed  and  SST  perturbations  have  a  strong  positive 
correlation  in  regions  of  large  SST  gradients,  which 
occur  near  meandering  ocean  currents  (see  review  by 
Small  et  al.  2008).  Contemporaneous  near-global  satel¬ 
lite  measurements  of  surface  vector  winds  and  SST  have 
shown  that  these  small-scale  features  in  the  surface  wind 
field  are  particularly  prevalent  near  large  midlatitude 
ocean  currents.  Previous  studies  have  only  investigated 
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Table  1.  Slopes  of  the  linear  relationships  between  the  wind  stress  curl  and  crosswind  SST  gradient  (ac)  and  the  wind  stress  divergence 
and  downwind  SST  gradient  (ao)  listed  in  previous  studies  over  various  regions  in  units  of  N  per  °C.  Acronyms  used  here  include: 
National  Center  for  Atmospheric  Research  (NCAR)  Community  Climate  System  Model  3.0  (CCSM3.0);  Model  for  Interdisciplinary 
Research  on  Climate  3.2  (MIROC3.2);  Scripps  Coupled  Ocean-Atmosphere  Regional  (SCOAR)  model;  and  Naval  Research  Laboratory 
Coupled  Ocean- Atmosphere  Mesoscale  Prediction  System  (CO AMPS)  model.  The  ratio  aJoiD  is  shown  in  the  rightmost  column. 


Study 

Region 

Source 

ac 

OLD 

olJold 

Chelton  et  al.  (2001) 

Tropical  Pacific 

3°N-1°S 

QuikSCAT-TMI 

0.0075 

0.0135 

0.56 

1°S-5°S 

0.0112 

0.0247 

0.45 

O’Neill  et  al.  (2003) 

Southern  Ocean 

QuikSCAT-Reynolds 

0.0068 

0.0124 

0.55 

Chelton  et  al.  (2004) 

Southern  Ocean 

QuikSCAT-AMSR 

0.0117 

0.0165 

0.71 

Tropical  Pacific 

0.0084 

0.0188 

0.45 

Kuroshio 

0.0057 

0.0096 

0.59 

Gulf  Stream 

0.0057 

0.0109 

0.52 

Chelton  (2005) 

Tropical  Pacific 

QuikSCAT-TMI 

0.0087 

0.0145 

0.60 

0.0080 

0.0142 

0.56 

O’Neill  et  al.  (2005) 

Agulhas 

QuikSCAT-AMSR 

0.0213 

0.0284 

0.75 

Maloney  and  Chelton  (2006) 

Agulhas 

QuikSCAT-AMSR 

0.0126 

0.0173 

0.73 

ECMWF 

0.0041 

0.0077 

0.53 

NCAR  CCSM3.0 

0.0041 

0.0088 

0.47 

MIROC3.2 

0.0074 

0.0131 

0.56 

Kuroshio 

QuikSCAT-AMSR 

0.0044 

0.0070 

0.63 

ECMWF 

0.0016 

0.0040 

0.40 

NCAR  CCSM3.0 

0.0032 

0.0065 

0.49 

MIROC3.2 

0.0021 

0.0059 

0.36 

Seo  et  al.  (2007) 

Tropical  Pacific 

Coastal  California 

SCOAR 

0.009 

0.015 

0.60 

36°-43°N 

0.002 

0.005 

0.40 

32°-39°N 

0.006 

0.008 

0.75 

Chelton  et  al.  (2007) 

California  current  system 

QuikSCAT-AMSR 

0.0213 

0.0301 

0.71 

Haack  et  al.  (2008) 

California  current  system 

QuikSCAT  CO  AMPS 

0.0182 

0.0260 

0.70 

0.0157 

0.0193 

0.81 

SST-induced  changes  in  wind  speed  across  SST  fronts 
and  not  wind  direction.  The  main  objectives  of  this  study 
are  (i)  to  determine  how  small-scale  SST  perturbations 
affect  both  surface  wind  speed  and  direction  from  6+ 
years  of  satellite  observations  of  vector  surface  winds 
and  SST  and  (ii)  to  determine  their  cumulative  effects 
on  the  surface  curl  and  divergence  fields  over  the  SST 
fronts  associated  with  four  midlatitude  ocean  currents: 
the  Kuroshio  Extension  in  the  North  Pacific,  the  Gulf 
Stream  in  the  North  Atlantic,  the  Brazil-Malvinas  Con¬ 
fluence  in  the  South  Atlantic,  and  the  Agulhas  Return 
Current  in  the  Indian  Ocean  sector  of  the  Southern 
Ocean. 

Spatial  variations  in  surface  wind  speed  associated 
with  small-scale  SST  variability  develop  surface  curl  and 
divergence  perturbations  with  magnitudes  comparable 
to  the  large-scale  curl  and  divergence  (Chelton  et  al. 
2004;  O’Neill  et  al.  2003,  2005;  Chelton  et  al.  2007). 
Satellite  observations  of  these  curl  and  divergence  per¬ 
turbations  are  remarkably  well  correlated  with  small- 
scale  perturbations  in  the  crosswind  and  downwind 
components  of  the  SST  gradient,  respectively  (Chelton 
et  al.  2001,  2004,  2007;  Chelton  2005;  O’Neill  et  al.  2003, 


2005).  The  surface  wind  stress  curl  and  divergence  fields 
are  related  linearly  to  the  crosswind  and  downwind  SST 
gradients,  respectively.  These  curl  and  divergence  de¬ 
pendencies  are  simulated  to  varying  degrees  in  numer¬ 
ical  weather  predication  and  climate  models  (Maloney 
and  Chelton  2006;  Haack  et  al.  2008)  and  in  regional 
mesoscale  models  (Seo  et  al.  2007;  Spall  2007b). 

In  all  past  studies,  the  curl  response  to  crosswind  SST 
gradients  has  been  found  to  be  significantly  weaker  than 
the  divergence  response  to  downwind  SST  gradients,  as 
measured  by  the  slopes  ac  and  ao,  respectively,  of  the 
linear  statistical  relations  (Table  1).  The  ratio  of  these 
slopes  aJaD  varies  between  about  0.4  and  0.75.  Addi¬ 
tionally,  satellite  observations  have  shown  strong  sea¬ 
sonal  variability  in  ac  and  a^  (O’Neill  et  al.  2005; 
Maloney  and  Chelton  2006;  O’Neill  et  al.  2009,  man¬ 
uscript  submitted  to  /.  Climate,  hereafter  OCE)  and 
strong  geographical  variability  between  ocean  basins 
(Chelton  et  al.  2004;  OCE)  and  within  individual  regions 
(O’Neill  et  al.  2005).  From  a  mesoscale  modeling  study 
over  the  Gulf  Stream,  Spall  (2007b)  noted  that  ac  has 
a  significant  quadratic  dependence  on  the  large-scale 
geostrophic  wind  speed.  This  is  shown  by  OCE  to  result 
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from  the  nonlinearity  between  SST-induced  surface  wind 
speed  and  wind  stress  perturbations. 

Beyond  what  is  known  about  the  dynamics  of  the 
near-surface  wind  response  to  small-scale  SST  pertur¬ 
bations,  as  discussed  in  section  2,  little  is  known  about 
the  dynamics  underlying  the  curl  and  divergence  re¬ 
sponses  to  SST  gradients.  The  differences  in  the  curl  and 
divergence  responses  reveal  an  opportunity  to  better 
understand  how  the  marine  atmospheric  boundary  layer 
(MABL)  responds  to  small-scale  SST  perturbations. 

In  this  study,  surface  wind  speed  and  direction  are 
obtained  from  the  SeaWinds  scatterometer  onboard  the 
Quick  Scatterometer  (QuikSCAT)  satellite  and  SST  is 
obtained  from  the  Advanced  Microwave  Scanning  Ra¬ 
diometer  for  Earth  Observing  System  (EOS)  (AMSR-E) 
Aqua  satellite,  as  discussed  in  section  3.  The  statistical 
responses  of  spatially  high-pass  filtered  surface  wind 
speed  and  direction  as  functions  of  SST  are  shown  in 
section  4.  A  relatively  simple  statistical  relation  is  found 
between  surface  wind  speed  and  SST  despite  the  com¬ 
plicated  force  balances  involved  in  the  SST-induced 
MABL  response  described  by  previous  studies.  In  sec¬ 
tion  5,  the  QuikSCAT  vorticity  and  divergence  of  the 
surface  wind  field  are  shown  to  depend  on  the  crosswind 
and  downwind  SST  gradients  in  a  manner  analogous  to 
the  dependencies  of  the  wind  stress  curl  and  divergence 
considered  in  our  previous  studies.  While  the  vorticity 
and  divergence  responses  to  SST  gradients  have  been 
interpreted  previously  in  terms  of  SST  effects  on  surface 
wind  stress  magnitude  (or  wind  speed),  we  show  that 
SST  also  modifies  the  surface  wind  direction,  thus  also 
contributing  to  SST-induced  changes  in  the  surface 
vorticity  and  divergence.  These  SST-induced  modifica¬ 
tions  of  wind  direction  are  shown  to  be  responsible  for 
the  differing  vorticity  and  divergence  responses  to  SST. 
We  close  in  section  6  with  a  discussion  of  the  implica¬ 
tions  of  the  results  of  this  study  for  coupled  mesoscale 
wind-SST  interactions. 

2.  Dynamical  explanations  for  mesoscale 
wind-SST  coupling 

Many  observational  and  modeling  studies  have  now 
been  undertaken  to  better  understand  the  MABL  re¬ 
sponse  to  SST  perturbations  over  various  regions  of  the 
World  Ocean.  The  response  of  surface  winds  to  meso¬ 
scale  SST  perturbations  is  the  culmination  of  adjustment 
processes  extending  throughout  the  depth  of  the  MABL 
and  appears  to  be  controlled  mainly  by  modification  of 
the  surface  heat  fluxes  as  air  flows  across  SST  fronts. 
Surface  heat  fluxes  are  enhanced  or  suppressed  depend¬ 
ing  on  whether  the  flow  is  from  cold  to  warm  water  or 
vice  versa.  Cross-frontal  surface  heating  perturbations 


produce  changes  in  the  MABL  hydrostatic  pressure 
gradients  and  to  the  vertical  turbulent  stress  divergence 
in  the  MABL  as  described  below,  ultimately  leading  to 
changes  in  near-surface  winds  coupled  to  the  mesoscale 
SST  perturbations. 

Across  SST  fronts,  variations  in  surface  heating  cause 
cross-frontal  pressure  gradients  near  the  surface  since 
a  cooler  MABL  over  cooler  water  forms  higher  surface 
pressures  relative  to  a  warmer  MABL  over  warmer 
water.  Thermally  induced  pressure  gradients  so  formed 
can  therefore  accelerate  near-surface  flow  across  SST 
isotherms  from  cooler  to  warmer  water  and  vice  versa. 
Using  an  analytical  model  of  the  boundary  layer, 
Lindzen  and  Nigam  (1987)  attributed  acceleration  of 
cross-equatorial  surface  flow  north  of  the  equatorial  cold 
tongue  in  the  eastern  tropical  Pacific  primarily  to  ther¬ 
mally  induced  pressure  gradients  associated  with  the  cold 
tongue  meridional  SST  gradient.  Even  though  some  of 
the  assumptions  of  their  analytical  model  are  not  gener¬ 
ally  valid  (e.g.,  Battisti  et  al.  1999;  Stevens  et  al.  2002; 
Small  et  al.  2003),  SST-induced  pressure  gradients  are 
indeed  strong  contributors  to  the  SST-induced  surface 
wind  response  (e.g.,  Wai  and  Stage  1989;  Warner  et  al. 
1990;  Small  et  al.  2003;  Cronin  et  al.  2003;  Mahrt  et  al. 
2004;  Bourras  et  al.  2004;  Song  et  al.  2004;  Small  et  al.  2005b; 
Song  et  al.  2006).  Using  mesoscale  numerical  simula¬ 
tions  over  the  equatorial  Pacific,  Small  et  al.  (2003) 
showed  that  near-surface  pressure  perturbations  form 
in  response  to  SST-induced  MABL  air  temperature  var¬ 
iations,  which  are  governed  thermodynamically  by  a  bal¬ 
ance  between  horizontal  temperature  advection  and 
surface  sensible  heat  fluxes.  Lrom  surface  pressure  obser¬ 
vations  over  the  equatorial  Pacific,  Cronin  et  al.  (2003) 
have  shown  that  SST-induced  surface  pressure  gradient 
perturbations  associated  with  the  northern  equatorial 
cold  tongue  are  of  a  magnitude  similar  to  those  gener¬ 
ated  by  mesoscale  model  simulations.  Hashizume  et  al. 
(2002)  have  also  suggested  that  deepening  and  shoaling 
of  the  MABL  associated  with  cross-frontal  turbulent 
mixing  variations  counteracts  this  thermodynamic  con¬ 
tribution  to  the  SST-induced  surface  pressure  variations 
in  what  has  been  referred  to  as  a  “backpressure”  effect. 
Finally,  while  MABL  pressure  gradients  contribute  to 
near-surface  flow  acceleration  perpendicular  to  SST 
isotherms,  it  is  not  clear  how  they  contribute  to  the 
generation  of  near-surface  vorticity  as  air  blows  parallel 
to  SST  isotherms,  a  feature  generally  observed  in  scat¬ 
terometer  wind  fields  globally. 

Wallace  et  al.  (1989)  and  Hayes  et  al.  (1989)  argued 
that  vertical  turbulent  mixing  of  momentum  from  aloft 
to  the  surface  was  more  consistent  with  the  changes  in 
vertical  wind  shear  and  the  acceleration  of  surface  winds 
occurring  across  the  northern  edge  of  the  equatorial 
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Pacific  cold  tongue  than  with  thermally  induced  pressure 
gradients.  Earlier  aircraft  observations  over  the  north 
wall  of  the  Gulf  Stream  by  Sweet  et  al.  (1981)  also  showed 
that  turbulent  mixing  of  momentum  from  aloft  to  the 
surface  also  played  a  significant  role  in  near-surface  wind 
speed  variations.  With  this  mechanism,  SST-induced  sur¬ 
face  heating  modifies  the  static  stability  of  the  boundary 
layer,  enhancing  the  vertical  turbulent  mixing  of  momen¬ 
tum  as  air  blows  from  cool  to  warm  water  and  reducing 
it  as  air  blows  from  warm  to  cool  water.  Observations 
of  the  near-surface  vertical  turbulent  momentum  flux 
throughout  the  world’s  oceans  have  generally  shown  con¬ 
siderable  variation  upon  crossing  SST  fronts  and  have 
supported  this  hypothesis  (Mey  and  Walker  1990;  Freihe 
et  al.  1991;  Bond  1992;  Jury  1994;  Rouault  and  Lutjeharms 
2000;  Mahrt  et  al.  2004;  de  Szoeke  et  al.  2005;  Tokinaga 
et  al.  2006). 

Modeling  studies  of  the  surface  wind  response  to  me- 
soscale  SST  perturbations  have  been  less  consistent 
about  the  role  that  the  vertical  turbulent  stress  divergence 
plays  in  the  SST-induced  surface  wind  response.  Some 
investigators  find  a  significant  role  (Wai  and  Stage  1989; 
de  Szoeke  and  Bretherton  2004;  Bourras  et  al.  2004; 
Song  et  al.  2004;  Skyllingstad  et  al.  2006;  Thum  2006; 
Song  et  al.  2009;  O’Neill  et  al.  2010)  while  others  found 
less  significant  roles  (Small  et  al.  2003;  Bourras  et  al. 
2004;  Small  et  al.  2005a,b;  Song  et  al.  2006;  Samelson 
et  al.  2006;  Spall  2007b).  Additionally,  the  contribution 
of  surface  ocean  currents  to  the  surface  stress  may  play 
a  small  but  nonnegligible  role  in  the  near-surface  verti¬ 
cal  turbulent  stress  divergence  (Song  et  al.  2006;  Liu 
et  al.  2007).  Idealized  large-eddy  simulations  of  the  flow 
across  sharp  SST  fronts  by  Skyllingstad  et  al.  (2006) 
showed  that  on  spatial  scales  of  1-20  km,  which  are 
somewhat  smaller  than  considered  here,  near-surface 
wind  speed  variations  are  influenced  predominantly  by 
SST-induced  cross-frontal  variations  in  the  vertical  tur¬ 
bulent  mixing  of  momentum,  consistent  with  the  Wallace 
et  al.  (1989)  and  Hayes  et  al.  (1989)  hypotheses.  Finally, 
weak  coupling  of  surface  winds  to  SST  in  the  European 
Centre  for  Medium-Range  Weather  Forecasts  (ECMWF) 
model  over  the  Agulhas  Return  Current  (Maloney  and 
Chelton  2006)  was  shown  by  Song  et  al.  (2009)  to  be 
attributable  primarily  to  a  weak  response  of  the  pa¬ 
rameterized  vertical  turbulent  stress  divergence  to  SST- 
induced  surface  heating  perturbations. 

Small  et  al.  (2005a)  envisioned  a  quite  different  role  of 
the  vertical  turbulent  stress  divergence  in  a  numerical 
simulation  of  the  near-surface  wind  response  over  the 
eastern  tropical  Pacific.  There,  it  was  found  that  pressure 
gradients  accelerated  the  near-surface  flow  across  the 
cold  tongue  SST  front,  and  after  some  distance,  the 
vertical  turbulent  stress  divergence  acted  as  a  frictional 


drag  on  this  pressure-driven  flow.  In  this  “pressure-drag” 
mechanism,  turbulence  acts  to  reduce  the  surface  winds 
in  an  opposite  manner  to  that  hypothesized  by  Wallace 
et  al.  (1989)  and  Hayes  et  al.  (1989).  This  pressure-drag 
mechanism  was  found  to  operate  in  a  numerical  simu¬ 
lation  over  the  Agulhas  Return  Current  (O’Neill  et  al. 
2010),  but  only  above  about  150  m  from  the  surface; 
below  this  level,  the  near-surface  winds  were  influenced  by 
the  turbulent  mixing  of  momentum  in  accordance  with  the 
Wallace  et  al.  (1989)  and  Hayes  et  al.  (1989)  mechanism. 

Samelson  et  al.  (2006)  argued  that  changes  in  boundary 
layer  depth  somewhat  downwind  of  the  sharpest  SST 
gradients,  rather  than  the  vertical  redistribution  of  mo¬ 
mentum,  could  act  to  change  the  surface  stress.  This 
mechanism  acts  through  a  balance  of  a  fixed  large-scale 
pressure  gradient  and  the  vertical  turbulent  stress  di¬ 
vergence  integrated  vertically  over  the  depth  of  the 
boundary  layer  such  that  the  ratio  of  surface  stress  to 
boundary  layer  height  tIH  is  equal  to  a  fixed  large-scale 
pressure  gradient  integrated  over  the  depth  of  the  bound¬ 
ary  layer.  This  balance  was  found  to  be  consistent  with 
the  cross-frontal  changes  in  surface  stress  and  boundary 
layer  height  away  from  the  steepest  SST  gradients  in 
the  two-dimensional  model  simulation  of  Spall  (2007b). 
This  balance  assumes  that  the  SST-induced  responses  of 
the  vertical  turbulent  mixing  of  momentum,  thermally 
induced  pressure  gradients,  and  horizontal  advection 
are  not  significant  terms  in  the  MABL  momentum  bud¬ 
get  in  regions  directly  over  the  sharpest  SST  fronts  (see 
detailed  discussion  in  Small  et  al.  2008). 

At  midlatitudes,  Coriolis  accelerations  have  been 
shown  to  be  important  to  the  response  of  the  surface 
winds  to  small-scale  SST  perturbations  (e.g.,  Wai  and 
Stage  1989;  Bourras  et  al.  2004;  Thum  2006;  Song  et  al. 
2006;  Spall  2007b;  O’Neill  et  al.  2010).  Over  a  single  SST 
front.  Spall  (2007b)  has  found  that  SST-induced  turbu¬ 
lent  mixing  perturbations  induce  Coriolis  accelerations 
analogous  to  an  inertially  driven  lee  wave.  O’Neill  et  al. 
(2010)  has  found  evidence  of  a  similar  mechanism  in 
a  simulation  over  a  part  of  the  Southern  Ocean.  The 
strong  large-scale  mean  wind  speeds  typically  found 
over  midlatitude  ocean  fronts  also  lead  to  significant 
horizontal  advective  accelerations  that  strongly  influ¬ 
ence  the  near-surface  momentum  budget  (Song  et  al. 
2006;  Thum  2006;  Spall  2007b;  O’Neill  et  al.  2010). 

Some  debate  remains  as  to  the  physical  mechanisms 
involved  in  the  coupling  between  surface  winds  and 
mesoscale  SST  perturbations.  It  is  becoming  clear,  how¬ 
ever,  that  both  pressure  gradients  and  the  turbulent 
mixing  of  momentum  are  involved  significantly  in  the 
response.  There  are  also  likely  significant  regional  dif¬ 
ferences  in  the  large-scale  atmospheric  forcing  that  may 
alter  the  force  balance  involved  in  the  surface  wind 
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response.  We  show  below  that  there  are  comparatively 
simple  statistical  relationships  between  the  surface  winds 
and  SST  despite  the  rather  complicated  details  of  these 
physical  mechanisms. 

3.  Description  of  satellite  observations  and  methods 

This  study  utilizes  monthly-averaged  surface  vector 
wind  observations  from  the  QuikSCAT  scatterometer 
and  SST  from  the  AMSR-E  over  the  75-month  period 
June  2002- August  2008,  coinciding  with  the  first  6+ 
years  of  the  AMSR-E  geophysical  observational  data 
record  that  began  1  June  2002.  Each  of  these  datasets  is 
summarized  in  this  section. 

a.  QuikSCAT  wind  fields 

In  nonprecipitating  conditions,  scatterometers  infer 
surface  vector  winds  over  water  from  microwave  radar 
measurements  of  small-scale  surface  roughness  caused 
by  surface  wind  stress.  For  lack  of  abundant  direct  sur¬ 
face  stress  measurements,  scatterometer  measurements 
of  microwave  radar  backscatter  are  calibrated  to  buoy 
anemometer  measurements  of  the  so-called  equivalent 
neutral  stability  wind  at  10  m,  which  is  the  10-m  wind 
that  would  be  uniquely  related  to  the  observed  surface 
wind  stress  if  the  atmosphere  were  neutrally  stratified 
(Liu  and  Tang  1996).  The  relation  between  the  equiva¬ 
lent  neutral  stability  wind  vector  u  and  the  surface  wind 
stress  vector  t  is  t  =  poC^Eu,  where  po  is  the  surface  air 
density,  V  is  the  magnitude  of  u,  and  is  the  neutral 
stability  drag  coefficient. 

The  computations  throughout  this  paper  using 
QuikSCAT  winds  are  based  on  10-m  equivalent  neutral 
stability  winds.  On  average,  anemometer  measurements 
of  the  10-m  wind  speed  are  about  0.2  m  s~^  lower  than 
the  corresponding  neutral-stability  wind  speed  at  10  m 
(Mears  et  al.  2001;  see  also  Fig.  16  of  Chelton  and 
Freilich  2005)  because  the  atmospheric  surface  layer  is, 
on  average,  slightly  unstable  over  all  oceans.  Through 
comparisons  with  high-quality  buoy  anemometer  mea¬ 
surements,  QuikSCAT  wind  speed  measurement  errors 
have  been  estimated  to  be  about  1.7  m  s“^  (Chelton  and 
Freilich  2005).  The  wind  direction  accuracy  improves 
with  increasing  wind  speed;  for  wind  speeds  greater  than 
about  5  m  s“^,  the  direction  accuracy  of  individual  wind 
measurements  is  better  than  15°  (see  Fig.  8  in  Chelton 
and  Freilich  2005).  Additionally,  there  is  no  evidence  of 
SST-dependent  measurement  biases  in  the  QuikSCAT 
wind  observations  (Chelton  et  al.  2001;  Ebuchi  et  al. 
2002). 

The  spatial  derivatives  of  wind  components,  wind 
speed,  and  wind  direction  analyzed  in  this  study  were 
computed  “in-swath”  before  interpolating  to  a  regular 


spatial  grid.  This  avoids  complications  arising  from  the 
estimation  of  spatial  derivatives  near  swath  edges  where 
centered  first-difference  estimates  of  derivatives  are  com¬ 
puted  using  measurements  from  neighboring  swaths  sep¬ 
arated  in  time  by  several  hours  or  more.  The  various  wind 
and  derivative  wind  fields  of  interest  in  this  study  were 
then  constructed  on  a  0.25°  latitude  by  0.25°  longitude 
spatial  grid  by  fitting  in-swath  measurements  to  a  qua¬ 
dratic  surface  using  locally  weighted  regression  (“loess” 
smoothing;  Cleveland  and  Devlin  1988)  with  a  half  span  of 
80  km.  Based  on  the  filter  transfer  function  of  the  loess 
smoother  (Schlax  et  al.  2001;  Chelton  and  Schlax  2003), 
the  resulting  gridded  fields  have  a  spatial  resolution 
analogous  to  that  of  approximately  50-km  block  averages. 
The  spatial  derivative  fields  were  found  at  each  grid  point 
directly  from  the  regression  coefficients  of  the  quadratic 
surface.  The  wind  fields  were  then  averaged  at  monthly 
intervals  from  the  gridded  swath  data. 

b.  AMSR-E  SST  fields 

The  detailed  results  of  this  study  would  not  be  possible 
without  the  near- all-weather  microwave  observations  of 
SST  over  middle-  and  high-latitude  regions  provided  by 
the  AMSR-E  since  June  2002  (Chelton  and  Wentz  2005). 
The  main  advantage  of  microwave  measurements  of  SST 
is  the  ability  to  measure  SST  through  nonprecipitating 
clouds,  which  are  essentially  transparent  to  microwave 
radiation.  Cloud  cover  generally  obscures  the  midlat¬ 
itude  ocean  regions  of  interest  in  this  study  more  than 
70%  of  the  time  (Chelton  and  Wentz  2005).  Large  re¬ 
gions  of  missing  data  often  occur  in  SST  fields  con¬ 
structed  from  infrared  satellite  measurements  over  the 
sharpest  ocean  fronts  and  eddies  where  strong  ocean- 
atmosphere  interactions  are  expected  to  occur  (see,  e.g.. 
Fig.  6  of  Chelton  and  Wentz  2005).  Indeed,  stratocu- 
mulus  and  cumulus  cloud  bands  often  form  over  sharp 
SST  fronts  as  a  consequence  of  some  of  the  same  MABL 
adjustment  processes  affecting  MABL  winds  of  inter¬ 
est  in  this  study  (e.g.,  Norris  and  lacobellis  2005).  While 
SST  measurements  by  the  Tropical  Rainfall  Measuring 
Mission  (TRMM)  Microwave  Imager  (TMI)  have  been 
available  since  December  1997  and  thus  encompass  the 
entire  QuikSCAT  mission,  they  are  only  made  equator- 
ward  of  40°  latitude,  leaving  most  of  the  midlatitude  re¬ 
gions  of  interest  in  this  analysis  unsampled. 

The  AMSR-E  measures  horizontally  and  vertically 
polarized  brightness  temperatures  at  six  microwave 
frequencies  across  a  1450-km  swath.  From  these  12  mi¬ 
crowave  brightness  temperatures,  SST  is  estimated  over 
most  of  the  global  oceans  using  a  physically  based  statis¬ 
tical  regression  algorithm  (Wentz  and  Meissner  2000).  The 
footprint  size  of  AMSR-E  measurements  of  SST  is  56  km 
and  the  rms  accuracy  of  individual  SST  measurements  is 
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about  0.4°C  (Chelton  and  Wentz  2005).  The  SST  fields 
were  gridded  onto  the  same  0.25°  grid  as  the  QuikSCAT 
surface  wind  fields  and  then  averaged  at  monthly  in¬ 
tervals. 

c.  Crosswind  and  downwind  SST  gradient 
computations 

As  discussed  in  the  introduction,  previous  studies  have 
shown  that  the  wind  stress  curl  and  divergence  fields 
are  found  to  be  related  linearly  to  the  crosswind  and 
downwind  components  of  the  SST  gradient,  respectively. 
Crosswind  and  downwind  gradients  in  this  study  were 
computed  within  each  QuikSCAT  measurement  swath 
using  curvilinear  natural  coordinates  (5,  n),  where  5-  and 
n  are  local  along-wind  and  crosswind  coordinates,  re¬ 
spectively  (e.g.,  Haltiner  and  Martin  1957;  Holton  1992). 
The  unit  vector  s  is  in  the  same  direction  as  u.  The 
positive  unit  vector  n  points  90°  counterclockwise  rela¬ 
tive  to  s.  The  crosswind  and  downwind  components  of 
the  SST  field  T{x,  y)  in  natural  coordinates  are  com¬ 
puted  from  Cartesian  derivatives  by 

dT  .  dT  ar  , 

—  =  -smil/ —  +  cosil/ —  and 
dn  dx  dy 

dT  dT  .  dT 

^  =  cosiIj^  +  sim/^^,  (1) 

ds  dx  dy 

where  (x,  y)  are  the  zonal  and  meridional  Cartesian 
coordinates,  respectively,  and  ip  is  the  counterclockwise 
surface  wind  direction  relative  to  the  x  axis. 

Crosswind  and  downwind  SST  gradients  were  computed 
here  on  a  swath-by-swath  basis  by  combining  wind  mea¬ 
surements  from  each  QuikSCAT  swath  with  SST  gradi¬ 
ents  computed  from  gridded  3-day  averaged  AMSR-E 
SST  fields  centered  on  each  day  of  the  QuikSCAT  mea¬ 
surement  swath.  Swaths  of  crosswind  and  downwind  SST 
gradients  were  then  averaged  at  monthly  intervals.  This 
procedure  minimizes  uncertainties  associated  with  com¬ 
puting  the  nonlinear  crosswind  and  downwind  SST  gra¬ 
dients  from  wind  and  SST  measurements  that  are  not 
collocated  spatially  and  temporally. 

d.  Spatial  high-pass  filtering 

SST-induced  mesoscale  wind  variability  was  isolated 
by  removing  large-scale  spatial  variability  using  a  mul¬ 
tidimensional  loess  smoothing  function  with  half-power 
filter  cutoffs  of  10°  latitude  X  20°  longitude  (Cleveland 
and  Devlin  1988).  This  filtering  effectively  removes  the 
large-scale  components  of  the  wind  and  SST  fields  that 
are  unrelated  to  the  wind-SST  coupling  of  interest  here. 
Wind  and  SST  fields  spatially  high-pass  filtered  in  this 
way  are  hereafter  referred  to  as  perturbation  fields.  To 
isolate  the  SST  influence  on  surface  winds  that  is  of  in¬ 


terest  here,  the  effects  of  synoptic  weather  variability 
were  mitigated  by  using  monthly-averaged  wind  and 
SST  fields.  It  is  noted  that  this  spatial  filtering  was  ap¬ 
plied  on  all  variables  after  computations  involving  the 
wind  and  SST  fields  were  performed,  including  the  non¬ 
linear  crosswind  and  downwind  gradient  computations. 
This  is  preferable  to  computing  the  crosswind  and  down¬ 
wind  gradient  quantities  from  monthly-averaged  fields.^ 

4.  Observations  of  SST-induced  surface  wind 
response 

a.  Surface  wind  speed  response  to  SST 

Maps  of  the  QuikSCAT  perturbation  surface  wind 
speed  and  the  AMSR-E  SST  fields  scalar-averaged  over 
the  75-month  analysis  period  are  shown  in  Fig.  1  for  each 
of  the  four  regions  considered  here.  Large  horizontal 
variations  in  surface  wind  speed  occur  near  meandering 
SST  fronts,  with  stronger  wind  speeds  over  warmer  water 
and  weaker  wind  speeds  over  cooler  water.  Typically, 
surface  wind  speed  variations  of  more  than  1-2  m  s~^ 
accompany  SST  variations  of  about  2°-4°C  over  cross- 
frontal  distances  of  0(100  km).  This  is  an  appreciable 
fraction  of  the  unfiltered  scalar-averaged  wind  speeds  in 
these  regions  (Fig.  2),  which  are  typically  8-13  m  s~^  in 
these  long-term  time  averages.  For  the  time-averaged 
perturbation  fields,  the  spatial  cross-correlation  between 
the  perturbation  wind  speed  and  the  perturbation  SST 
varies  between  0.7  over  the  Kuroshio  to  0.85  over  the 
South  Atlantic  (Table  2).  In  addition  to  the  quasi-stationary 
coupled  patterns  in  the  perturbation  wind  and  SST  fields 
apparent  in  these  long-term  time-averaged  maps,  there 
is  also  significant  coupled  variability  associated  with 
transient  SST  features  evident  from  satellite  data  over  the 
Kuroshio  (e.g.,  Nonaka  and  Xie  2003),  the  Gulf  Stream 
(e.g..  Park  and  Cornillon  2002;  Park  et  al.  2006),  and  all 
four  midlatitude  regions  considered  here  (White  and 
Annis  2003).  The  wind-SST  coupling  in  the  Kuroshio  re¬ 
gion  is  much  more  apparent  during  the  boreal  wintertime 


^  Crosswind  and  downwind  gradients  of  SST  in  our  previous 
studies  (Chelton  et  al.  2001,  2004;  O’Neill  et  al.  2003,  2005)  were 
computed  from  gridded,  vector-averaged  wind  stress  and  SST 
gradient  components  rather  than  being  computed  on  a  swath-by¬ 
swath  basis  and  then  averaging,  which  is  more  accurate  because  of 
the  nonlinearity  of  these  quantities.  As  a  result,  the  SST-induced 
wind  response  in  those  studies  was  somewhat  underestimated  from 
those  obtained  from  the  more  accurately  computed  quantities  used 
here.  Additionally,  the  definition  of  the  crosswind  SST  gradient 
used  in  those  studies  (VT  X  u)  •  ^  is  equal  to  -dTIdn  and  the 
downwind  SST  gradient  \T  •  u  is  equal  to  dT/ds,  where  u  is  a  unit 
vector  in  the  direction  of  the  surface  wind  vector  and  ^  is  a  unit 
vector  in  the  vertical  direction. 
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a)  Kuroshio 
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Spatially  High-Pass  Filtered  Surface  Wind  Speed  (m/s) 

Fig.  1.  Maps  of  the  perturbation  QuikSCAT  wind  speed  (colors)  and  AMSR-E  SST  (contours)  scalar-averaged 
over  the  period  June  2002-August  2008  for  the  four  regions  considered  in  this  study.  The  contour  interval  for  the 
perturbation  SST  is  0.5°C,  and  the  zero  contour  has  been  omitted  for  clarity.  Solid  and  dashed  contours  correspond  to 
positive  and  negative  values,  respectively.  The  spatial  high-pass  filter  removes  spatial  variability  with  wavelengths 
longer  than  20°  longitude  X  10°  latitude,  as  discussed  in  the  text. 


(e.g.,  Nonaka  and  Xie  2003;  OCE)  and  is  associated  with 
strong  seasonality  of  the  SST  fronts  in  this  region  (e.g., 
Nonaka  and  Xie  2003). 

The  mesoscale  response  of  the  wind  speed  to  SST  is 
quantified  statistically  by  bin-averaging  monthly  averages 
of  the  perturbation  QuikSCAT  wind  speed  as  a  function 
of  the  perturbation  AMSR-E  SST  over  the  75-month 
analysis  period  (Fig.  3).  Wind  speed  perturbations  are 
related  linearly  to,  and  correlated  positively  with,  me¬ 
soscale  SST  perturbations  over  all  four  regions.  The 
relationship  between  the  monthly-averaged  spatially  high- 
pass  filtered  wind  speed  V'  and  SST  T'  can  thus  be  ex¬ 
pressed  empirically  as 

V'=aJ',  (2) 

where  is  the  least  squares  estimate  of  the  slope  of 
these  linear  relations  and  represents  the  change  in  wind 
speed  per  unit  change  in  SST  (i.e.,  =  dV'IdT)  and  the 

primes  hereafter  represent  spatially  high-pass  filtered 
fields.  Notably,  varies  by  nearly  a  factor  of  2  geo¬ 
graphically,  from  0.27  m  s~^  per  °C  over  the  Gulf  Stream 
to  0.49  m  s~^  per  °C  over  the  Agulhas  Return  Current 


(see  Table  2).  Additionally,  is  much  larger  in  the  two 
Southern  Hemisphere  regions  compared  to  those  in  the 
Northern  Hemisphere.  The  larger  error  bars  at  large 
perturbation  SST  magnitudes  in  Fig.  3  are  due  to  the 
smaller  numbers  of  samples  in  these  bins.  This  is  par¬ 
ticularly  acute  over  the  Kuroshio,  where  the  monthly- 
averaged  SST  perturbations  tend,  on  average,  to  be 
smaller  in  magnitude  than  in  the  other  three  regions. 
The  geographical  differences  in  between  the  four 
regions  are  likely  due  to  geographic  differences  in  the 
vertical  structure  of  the  boundary  layer  and  large-scale 
forcing.  Seasonality  of  is  investigated  from  5  +  years 
of  satellite  data  over  these  four  regions  (and  the  eastern 
tropical  Pacific)  in  OCE. 

Finally,  we  note  that  for  SST  perturbations  warmer 
than  about  1.5°,  there  is  a  subtle  flattening  of  the  bin- 
averaged  V'  relative  to  the  linear  fit,  particularly  over 
the  South  Atlantic  and  Agulhas  Return  Current  regions. 
Interestingly,  this  flattening  does  not  occur  over  cool 
SST  perturbations.  This  difference  between  the  bin- 
averaged  V'  and  the  linear  fit  is  less  than  about  0.25  m  s“^, 
and  there  are  very  few  data  points  in  these  outer  bins,  as 
shown  by  the  histograms  in  Fig.  3.  Because  this  flattening 
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Fig.  2.  Maps  of  the  unfiltered  QuikSCAT  wind  speed  (colors)  and  AMSR-E  SST  (contours)  scalar-averaged  over 
the  period  June  2002-August  2008.  The  SST  contour  interval  is  2°C. 


of  V'  is  relatively  small  and  affects  only  a  very  small 
percentage  of  the  observations,  it  does  not  significantly 
affect  the  analysis  presented  here.  This  observation  is 
currently  the  subject  of  ongoing  analysis. 

b.  Wind  direction  response  to  SST 

The  wind  direction  over  these  four  regions  is  shown 
in  Fig.  4  by  the  streamlines  of  the  vector-averaged  sur¬ 
face  wind  from  QuikSCAT  for  the  75-month  period  con¬ 
sidered  here.  The  surface  flow  in  all  four  regions  has  a 
strong  westerly  component,  characteristic  of  typical  mid¬ 
latitude  surface  flow.  The  flow  immediately  downstream 
of  the  Asian  and  North  American  continents  has  a  sig¬ 
nificant  northerly  component  before  gradually  turning 
northward  far  offshore.  Over  the  Southern  Hemisphere 
regions,  the  flow  is  westerly  poleward  of  the  large-scale 
subtropical  anticyclones  located  near  30°-35°S. 

For  quantitative  purposes,  it  is  desirable  here  to  consider 
wind  direction  perturbations  analogous  to  the  spatially 
high-pass  filtered  wind  speed  perturbations  investigated 
in  section  4a.  However,  scalar  wind  direction  cannot  be 
spatially  high-pass  filtered  in  the  same  manner  as  wind 
speed  in  regions  straddling  absolute  jumps  in  wind  di¬ 
rection  exceeding  360°  and  where  the  wind  speed  ap¬ 


proaches  zero,  in  which  case  the  wind  direction  becomes 
undefined.  Both  cases  occur  intermittently  in  each  of  the 
four  regions  under  consideration.  To  avoid  these  diffi¬ 
culties,  mesoscale  wind  direction  perturbations  are  de¬ 
fined  here  as  the  counterclockwise  angle  between  the 
unfiltered  u  and  the  spatially  low-pass  filtered  u,  where 
the  tilde  represents  spatial  low-pass  filtering.  These  wind 
direction  perturbations  were  computed  from  vector- 
averaged  wind  components  at  monthly  intervals  only  at 
points  where  the  magnitude  of  the  monthly  vector- 
averaged  wind  speed  exceeded  1  m  s~^. 

Spatial  maps  of  the  so-defined  mesoscale  wind  direction 
perturbations  scalar-averaged  over  the  75-month  anal¬ 
ysis  period  are  shown  in  Fig.  5  along  with  contours  of  the 


Table  2.  The  (left)  coupling  coefficient  and  (right)  spatial 
correlation  coefficient  computed  from  the  monthly-averaged  per¬ 
turbation  QuikSCAT  wind  speed  and  AMSR-E  SST  fields. 


[m  s  ^  (°C) 

Correlation  coefficient 

Kuroshio 

0.30 

0.70 

Gulf  Stream 

0.27 

0.79 

South  Atlantic 

0.43 

0.87 

Agulhas 

0.49 

0.83 
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Fig.  3.  Binned  scatterplots  of  the  spatially  high-pass  filtered  QuikSCAT  wind  speed  {V')  as  a  function  of  the 
spatially  high-pass  filtered  AMSR-E  SST  (T).  The  bin  averages  were  computed  from  monthly-averaged  QuikSCAT 
wind  and  AMSR-E  SST  fields  over  the  75-month  analysis  period  using  all  available  data  points  in  the  regions  en¬ 
closed  in  Figs.  1  and  2.  The  points  and  error  bars  represent  the  means  and  standard  deviations  within  each  bin, 
respectively,  and  the  dashed  lines  are  linear  least  squares  fits  to  the  means  within  each  bin.  Below  each  scatterplot, 
a  histogram  of  the  spatially  high-pass  filtered  AMSR-E  SST  is  included  for  reference. 


AMSR-E  perturbation  SST.  Perhaps  the  most  striking 
feature  of  these  wind  direction  perturbations  is  the  pro¬ 
nounced  meridional  displacement  of  the  wind  direction 
perturbations  relative  to  those  of  SST.  In  the  Northern 
Hemisphere,  negative  wind  direction  perturbations  are 
shifted  north  of  warm  SST  perturbations  and  positive 
wind  direction  perturbations  are  shifted  north  of  cool  SST 
perturbations.  Conversely,  in  the  Southern  Hemisphere, 
positive  wind  direction  perturbations  are  shifted  south 
of  warm  SST  perturbations  and  negative  wind  direction 
perturbations  are  shifted  south  of  cool  SST  perturbations. 

Meridional  displacements  of  the  wind  direction  per¬ 
turbations  relative  to  those  of  SST  that  are  qualitatively 
evident  in  Fig.  5  preclude  the  use  of  simple  binned  av¬ 
erages  to  statistically  quantify  the  relationship  between 
wind  direction  and  SST  such  as  was  performed  with  the 
wind  speed  perturbations  in  section  4a  (Fig.  3).  Statis¬ 
tical  relationships  between  the  wind  direction  and  SST 


perturbations  can  nonetheless  be  obtained  here  through 
cross-spectral  analysis  of  the  wind  direction  and  SST 
perturbations  in  the  meridional  wavenumber  domain. 
Estimates  of  cross-spectral  statistics  (squared  coherence, 
transfer  function,  and  phase  spectra  as  functions  of  me¬ 
ridional  wavenumber)  were  computed  at  monthly  in¬ 
tervals  along  continuous  meridional  segments  within  each 
of  the  four  regions  and  then  ensemble  averaged  over  all 
75  months  and  segments  (Fig.  6). 

The  squared  coherence  is  largest  over  the  South  At¬ 
lantic  and  Agulhas  Return  Current,  where  perturbation 
SST  variability  accounts  for  10%-15%  of  the  meridional 
variance  in  the  perturbation  wind  direction.  Over  the 
Gulf  Stream  and  Kuroshio,  the  squared  coherence  is 
below  0.05.  Wind  direction  perturbations  are  therefore 
less  well  correlated  with  SST  compared  with  wind  speed. 
The  transfer  function  represents  an  estimate  of  the  cou¬ 
pling  coefficient  between  the  perturbation  wind  direction 
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Fig.  4.  Maps  of  the  unfiltered  AMSR-E  SST  (colors)  scalar-averaged  over  the  period  June  2002-August  2008 
for  the  four  regions  considered  in  this  study.  Overlaid  are  select  streamlines  of  the  surface  flow  computed  from  the 
75-month  vector-averaged  QuikSCAT  winds,  where  the  arrows  indicate  flow  direction. 


and  SST  as  a  function  of  meridional  wavenumber.  From 
the  transfer  function  in  Fig.  6,  the  linear  responses  of  the 
wind  direction  perturbations  are  roughly  2°  per  degree 
Celsius  change  in  perturbation  SST  over  all  four  re¬ 
gions  and  vary  little  with  meridional  wavenumber.  SST- 
induced  wind  direction  perturbations  are  thus  about  4°-8° 
for  a  2°^°C  change  in  perturbation  SST,  which  is  con¬ 
sistent  with  the  magnitude  of  the  wind  direction  pertur¬ 
bations  evident  in  Fig.  5.  The  meridional  phase  spectrum 
is  a  measure  of  the  meridional  phase  difference  as  a 
function  of  meridional  wavenumber.  Over  all  four  re¬ 
gions,  the  phase  spectra  vary  approximately  linearly  with 
wavenumber,  indicative  of  a  constant  meridional  dis¬ 
placement  between  the  wind  direction  and  SST  pertur¬ 
bations.  The  meridional  displacement  Lq  of  the  wind 
direction  perturbations  relative  to  the  SST  perturbations 
estimated  from  the  slope  of  the  phase  spectra  in  Fig.  6  is 
poleward  about  1°  latitude  over  all  four  regions,  consis¬ 
tent  with  the  spatial  lags  evident  visually  in  Fig.  5. 

The  influence  of  SST  on  wind  direction  has  been  shown 
qualitatively  from  previous  case  studies  over  the  Gulf 
Stream  using  mesoscale  atmospheric  models  (Song  et  al. 
2006;  Spall  2007b)  and  scatterometer  observations  of 
Gulf  Stream  rings  (Park  et  al.  2006)  and  over  the  Southern 


Ocean  from  model  simulations  (Thum  2006;  O’Neill  et  al. 
2010)  and  aircraft  observations  (e.g..  Jury  1994).  The  results 
in  these  studies  are  consistent  with  the  results  shown  here 
based  on  75  months  of  QuikSCAT  wind  and  AMSR-E  SST 
observations. 

These  observations  of  the  SST-induced  wind  direction 
perturbations  show  that  1)  wind  direction  perturbations 
are  shifted  poleward  by  about  1°  latitude  relative  to  the 
SST  perturbations  and  2)  in  the  Northern  (Southern) 
Hemisphere,  warm  SST  perturbations  induce  negative 
(positive)  wind  direction  perturbations  and  cool  SST  per¬ 
turbations  induce  positive  (negative)  wind  direction 
perturbations;  the  magnitude  of  these  wind  direction 
perturbations  is  about  2°  per  degree  Celsius  change 
in  perturbation  SST  over  all  four  regions. 

The  wind  vector  responses  to  meanders  in  SST  fronts 
are  summarized  schematically  in  Fig.  7  for  each  hemi¬ 
sphere.  For  illustrative  purposes,  the  meridional  dis¬ 
placement  of  the  wind  direction  perturbations  relative 
to  the  SST  front  are  not  shown  here  because  they  are 
small  relative  to  the  overall  meridional  extent  of  the  SST 
perturbations,  which  is  typically  3°-5°  latitude.  Note  that 
for  surface  flow  from  cool  to  warm  SST,  the  surface 
winds  turn  anticyclonically,  whereas  for  surface  flow 
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Fig.  5.  As  in  Fig.  1,  but  for  the  QuikSCAT  perturbation  wind  direction  scalar  averaged  over  the  period  June  2002- 
August  2008.  As  discussed  in  the  text,  wind  direction  perturbations  are  defined  here  as  the  counterclockwise  angle 
between  the  unfiltered  and  spatially  low-pass  filtered  surface  wind  vectors. 


from  warm  to  cool  SST,  the  surface  winds  turn  cycloni- 
cally.^  Although  the  coupling  between  wind  direction 
and  SST  is  weaker  than  the  coupling  between  wind 
speed  and  SST,  it  nonetheless  has  important  conse¬ 
quences  for  the  surface  vorticity  and  divergence  fields, 
as  will  be  discussed  below.  Small  changes  in  wind  direction 
over  short  distances  produce  relatively  large  vorticity  and 
divergence  perturbations,  as  demonstrated  below. 


5.  SST-induced  variability  in  the  surface  vorticity 
and  divergence  fields 


The  separate  effects  of  spatial  gradients  in  wind  speed 
and  direction  on  the  vorticity  and  divergence  fields  are 
ascertained  by  expressing  the  vorticity  and  divergence  in 
curvilinear  natural  coordinates  (e.g.,  Haltiner  and  Martin 
1957): 


V 


(3) 


V  u'  = 


(4) 


where  (s,  n)  are  the  local  downwind  and  crosswind  co¬ 
ordinates,  respectively.  Equation  (3)  partitions  the  per¬ 
turbation  vorticity  field^  into  the  difference  between 
the  local  crosswind  speed  gradient  and  a  term  related 
to  the  local  downwind  gradient  in  flow  direction.  Likewise, 
Eq.  (4)  partitions  the  perturbation  divergence  field  into  the 
sum  of  the  downwind  speed  gradient  and  a  term  related 
to  the  crosswind  gradient  in  flow  direction.  The  direction 
gradient  term  in  the  vorticity  is  related  to  the  radius  of 
curvature  of  surface  streamlines,  whereas  the  direction 
gradient  term  in  the  divergence  is  related  to  spread¬ 
ing  or  contracting  of  streamlines  in  the  crosswind  direc¬ 
tion,  which  is  often  referred  to  as  the  diffluent  tendency 
(Haltiner  and  Martin  1957).  Decomposing  the  vorticity 
and  divergence  in  this  manner  allows  investigation 
of  the  separate  effects  of  lateral  shear  and  downwind 


^  Anticy clonic  rotation  is  clockwise  in  the  Northern  Hemisphere 
and  counterclockwise  in  the  Southern  Hemisphere.  Cyclonic  ro¬ 
tation  is  the  opposite  of  this. 


^  The  horizontal  vorticity  should  be  expressed  formally  as 
(V  X  u)  •  k.  The  expression  V  X  u  is  used  as  shorthand  for  this  to 
simplify  the  mathematical  notation  throughout  the  rest  of  this 
analysis. 
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Fig.  6.  Cross-spectral  statistics  of  the  perturbation  wind  direction  dependence  on  the  perturbation  SST  as 
a  function  of  meridional  wavenumber:  (top)  squared  coherence,  (middle)  transfer  function,  and  (bottom)  phase 
spectra.  Only  cross-spectral  estimates  at  wavenumbers  with  squared  coherences  above  the  95%  confidence  level  are 
shown.  The  dashed  lines  in  the  phase  spectra  plots  are  linear  least  squares  fits  of  the  phase  estimates  with  slopes  as 
indicated.  These  are  used  to  compute  estimates  of  the  meridional  displacement  Lq  between  wind  direction  and  SST 
perturbations,  which  are  indicated  in  each  panel;  the  units  of  Lq  shown  are  in  degrees  latitude. 


rotation  on  the  vorticity  and  the  separate  effects  of 
downwind  acceleration  and  crosswind  rotation  on  the 
divergence. 

Crosswind  and  downwind  wind  speed  gradients  were 
computed  in-swath  using 


'  .  dV  ,  dV\ 
—smi// —  +  cosi/j  — 

^  dx  dy  J 


'  ,  dV  .  dV\ 

cosi/j —  +  Sint//  — 
dx  dyj 


and 
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Fig.  7.  Summary  schematic  of  the  vector  wind  response  to  meanders  along  an  extratropical  SST  front,  represented 
here  by  the  solid  curves,  as  deduced  from  the  wind  speed  and  direction  dependencies  on  SST.  Whereas  the  wind 
speed  response  to  SST  (as  represented  by  the  relative  length  of  the  vectors)  is  the  same  for  both  hemispheres,  the 
wind  direction  response  to  SST  (as  represented  by  the  relative  turning  of  the  vectors)  differs  in  sign  between  the 
Northern  and  Southern  Hemispheres. 


The  crosswind  and  downwind  direction  gradients  were 
computed  from  the  gridded  vorticity,  divergence,  and 
speed  gradient  fields  using 


dn^ 


ds 


and 


V  X  u  +  —  . 

dn 


a.  Spatial  gradients  of  perturbation  wind  speed 
and  SST 

From  the  linear  statistical  relationship  between  V  and 
T  expressed  by  Eq.  (2),  it  follows  that  crosswind  and 
downwind  gradients  of  wind  speed  should  depend,  re¬ 
spectively,  on  the  crosswind  and  downwind  components 
of  the  SST  gradient  field  such  that 


and 


(5) 


(6) 


where  and  are  the  coupling  coefficients  for  the 
speed  gradient  responses  to  the  SST  gradients.  The 
QuikSCAT  crosswind  and  downwind  speed  gradient 
components  are  indeed  related  linearly  to  the  crosswind 
and  downwind  SST  gradient  components,  respectively, 
for  the  75-month  period  of  interest  here  (Fig.  8).  More¬ 


over,  within  each  region  and  are  nearly  equal  to 
from  Fig.  3.  The  scalar  wind  speed  and  speed  gradi¬ 
ent  responses  to  mesoscale  SST  perturbations  therefore 
appear  to  be  manifestations  of  the  same  physical  mech¬ 
anisms.  This  strong  coupling  between  the  crosswind 
and  downwind  speed  and  SST  gradients  is  further  illus¬ 
trated  in  the  maps  of  the  perturbation  crosswind  and 
downwind  speed  and  SST  gradients  averaged  over  the 
75-month  period  (Fig.  9).  The  corresponding  cross¬ 
correlation  coefficients  between  these  time-averaged 
fields  are  shown  in  Table  3;  they  are  comparable  to  those 
of  the  wind  speed  and  SST  perturbations  presented  in 
Table  2  for  the  Southern  Hemisphere,  while  being  some¬ 
what  weaker  in  the  Northern  Hemisphere.  Note  that  in 
the  core  of  the  Gulf  Stream  region  near  37°N,  70°W, 
there  is  little  signal  in  the  downwind  speed  gradients 
even  though  the  downwind  SST  gradients  are  strong. 
This  may  be  due  to  the  effects  of  the  strong  Gulf  Stream 
ocean  surface  currents  on  the  scatterometer  winds,  as 
argued  in  Fig.  5  of  Chelton  et  al.  (2004). 

The  linear  relationships  between  the  perturbation  wind 
speed  and  SST  fields  can  be  investigated  further  by  ex¬ 
pressing  the  perturbation  crosswind  and  downwind  com¬ 
ponents  of  the  SST  gradient  as 

=  Mj.smO'  and 
=  Mj.cos6',  (7) 
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Fig.  8.  As  in  Fig.  3,  but  for  the  QuikSCAT  perturbation  crosswind  speed  gradient  (dVIdn)'  binned  as  a  function  of 
the  crosswind  SST  gradient  {dTIdn)'  (black)  and  the  QuikSCAT  perturbation  downwind  speed  gradient  {dVIds)' 
binned  as  a  function  of  the  downwind  SST  gradient  {dTIds)'  (gray). 


where  the  perturbation  SST  gradient  magnitude  Mr  and 
the  angle  6'  are  defined  such  that 


M^  = 


and 


6'  =  tan  ^ 


{dT/dn) 
{dTIds)  ' 


(8) 


Geometrically,  6'  closely  approximates  the  angle  be¬ 
tween  the  surface  streamlines  of  the  unfiltered  surface 
winds  and  the  perturbation  SST  gradient  vector  \T 
(justification  of  this  physical  interpretation  of  is  given 
in  the  appendix).  From  consideration  of  the  simple  lin¬ 
ear  relationship  between  the  perturbation  wind  speed 
and  SST  in  Eq.  (2),  the  largest  downwind  speed  gradi¬ 
ents  should  occur  as  surface  winds  blow  across  pertur¬ 
bation  SST  isotherms  such  that  the  wind  speed  increases 
as  winds  blow  from  cool  to  warm  water  when  9'  =  0°, 
resulting  in  (dU/d^)'  >  0  in  accord  with  Eq.  (6),  and 
decreases  as  winds  blow  from  warm  to  cool  water  when 
6'  =  ±180°,  resulting  in  (dV/ds)'  <  0.  Likewise,  the 
largest  crosswind  speed  gradients  should  occur  as  the 


surface  winds  blow  along  perturbation  SST  isotherms 
(when  9'  =  ±90°). 

The  angular  dependencies  of  the  perturbation  crosswind 
and  downwind  speed  gradients  on  9'  are  shown  in  Fig.  10. 
The  crosswind  and  downwind  speed  gradients  closely 
follow  sine  and  cosine  functions,  respectively,  of  the  angle 
9',  in  agreement  with  the  hypothesized  dependencies 
relating  the  wind  speed  gradients  to  the  orientation  of 
the  surface  streamlines  relative  to  the  perturbation  SST 
isotherms.  Close  inspection  of  Fig.  10  reveals  small  but 
consistent  phase  shifts  from  the  expected  pure  cosine 
curves.  From  these  statistical  results,  we  may  thus  express 
the  crosswind  and  downwind  speed  gradients  as 

=  sin(6i'  +  and  (9) 

(10) 

where  and  are  the  slopes  of  the  speed  gradient 
responses  to  the  SST  gradients  shown  in  Fig.  8  and 
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crosswind  SST  gradient,  (bottom  maps  in  each  panel)  The  perturbation  downwind  speed  gradient  (colors)  overlaid 
with  contours  of  the  downwind  SST  gradient.  The  data  in  each  map  were  averaged  over  the  75-month  period  June 
2002-August  2008.  The  contour  interval  for  the  perturbation  crosswind  and  downwind  SST  gradients  is  0.2°C  per 
100  km,  and  the  zero  contour  has  been  omitted  for  clarity.  Solid  and  dashed  contours  correspond  to  positive  and 
negative  values  of  the  SST  gradient,  respectively. 


and  (/)  are  the  phase  shifts  of  the  sinusoidal  responses 
of  the  speed  gradients  to  6'  shown  in  Fig.  10  (positive  for 
counterclockwise  rotation).  The  values  of  the  ampli¬ 
tudes  and  phase  shifts  are  shown  in  Table  4.  Within  each 
region,  and  are  nearly  equal.  Note  that  (/)^p^  is 
essentially  0°  over  all  four  regions  whereas  is  8°-17° 
with  opposing  signs  between  hemispheres. 

b.  SST-induced  speed  and  direction  gradient  effects 
on  vorticity  and  divergence 

The  perturbation  vorticity  and  divergence  fields  binned 
as  functions  of  the  perturbation  crosswind  and  downwind 
SST  gradients,  respectively,  are  shown  in  Fig.  11.  The 
perturbation  vorticity  and  divergence  fields  are  highly 
correlated  with  small-scale  perturbations  in  the  SST 
gradient  field,  as  demonstrated  in  our  previous  studies 
using  the  wind  stress  curl  and  divergence  fields  (Chelton 
et  al.  2001,  2004,  2007;  O’Neill  et  al.  2003,  2005).  To  em¬ 
phasize  the  difference  between  the  vorticity  and  diver¬ 
gence  responses,  and  for  consistency  with  our  previous 


studies,  the  vorticity  is  binned  as  a  function  of  the  neg¬ 
ative  crosswind  SST  gradient  —(BTIdn)'  in  Fig.  11.  In 
agreement  with  these  previous  studies  shown  in  Table  1, 
the  magnitude  of  the  slope  of  the  divergence  response  to 
downwind  SST  gradients  is  significantly  larger  than  that 
of  the  vorticity  response  to  crosswind  SST  gradients. 
Additionally,  the  variability  within  each  bin,  as  repre¬ 
sented  by  the  standard  deviation  within  each  bin,  is  much 
smaller  for  the  speed  gradient  components  in  Fig.  8  than 
for  the  vorticity  and  divergence  in  Fig.  11.  This  is  evi¬ 
dently  because  contributions  of  the  SST-induced  wind 
direction  responses  to  the  vorticity  and  divergence  are 
noisier  than  the  wind  speed  responses. 

Given  the  results  of  section  3a  and  Eqs.  (3)  and  (4),  if 
the  surface  vorticity  and  divergence  were  dependent 
only  on  the  crosswind  and  downwind  speed  gradients, 
respectively,  then  the  slopes  of  the  straight  line  fits  in 
Fig.  11  should  be  equal  in  magnitude  but  opposite  in  sign. 
Evidently,  the  direction  gradient  terms  in  Eqs.  (3)  and 
(4)  reduce  the  overall  vorticity  responses  to  crosswind 
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Table  3.  Cross-correlation  coefficients  between  the  various 
perturbation  derivative  wind  fields  and  the  perturbation  crosswind 
and  downwind  SST  gradient  components  as  indicated.  These  cor¬ 
relation  coefficients  were  computed  from  the  derivative  wind  and 
SST  fields  averaged  over  the  75-month  analysis  period  at  monthly 
intervals. 
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SST  gradients  and  enhance  the  divergence  responses  to 
downwind  SST  gradients.  This  point  is  addressed  in  de¬ 
tail  in  section  4c. 

Besides  these  differences,  there  are  significant  differ¬ 
ences  in  the  vorticity  and  divergence  dependencies  on  O' 
(Fig.  12)  from  that  expected  solely  from  consideration  of 
SST-induced  wind  speed  perturbations.  In  our  previous 
work,  we  hypothesized  that  the  surface  curl  should  de¬ 
pend  on  the  sine  of  O'  whereas  the  divergence  should 
depend  on  the  cosine  of  O'.  Inspection  of  the  vorticity 
and  divergence  dependencies  on  in  Fig.  12,  however, 
reveals  that  the  vorticity  and  divergence  do  not  depend 
exactly  on  the  sine  and  cosine,  respectively,  of  O'  but 
rather  are  phase  shifted  (Fig.  12).  The  values  of  the 
phase  shifts  are  labeled  on  Fig.  12.  These  phase  shifts 
have  absolute  magnitudes  of  approximately  25°  for  the 
vorticity  and  10°  for  the  divergence  in  all  four  regions 
but  with  opposite  signs  in  each  hemisphere.  As  addressed 
below,  these  phase  shifts  and  the  amplitude  differences  in 
Fig.  11  (also  evident  from  the  amplitudes  of  the  sinusoids 
in  Fig.  12)  can  be  related  directly  to  the  influence  of  SST 
on  the  direction  gradient  terms  in  Eqs.  (3)  and  (4). 

The  dependence  of  the  direction  gradient  terms  on 
SST  can  be  ascertained  conceptually  from  the  schematic 
in  Fig.  7.  When  the  surface  flow  is  from  cool  to  warm 
SST,  the  surface  winds  turn  anticyclonically,  whereas  for 
surface  flow  from  warm  to  cool  SST  the  surface  winds 
turn  cyclonically.  This  observation  suggests  that  the 
downwind  gradient  in  wind  direction  V(diplds)  should 
depend  on  the  downwind  SST  gradient  and  hence  on  the 
cosine  of  O'.  Additionally,  crosswind  direction  gradients 
V(diljldn)  should  form  as  the  wind  blows  along  pertur¬ 


bation  SST  isotherms  in  association  with  the  crosswind 
SST  gradient  and  will  thus  depend  on  the  sine  of 

The  time-averaged  maps  of  the  perturbation  downwind 
direction  and  SST  gradients  and  the  perturbation  cross- 
wind  direction  and  SST  gradients  in  Fig.  13  indeed  show 
a  correlation  between  the  direction  gradient  and  SST 
gradients  consistent  with  these  hypothesized  dependen¬ 
cies.  The  poleward  displacement  of  the  wind  direction 
perturbations  relative  to  those  of  SST  discussed  in  section 
4b  is  also  evident  in  these  maps.  The  correlations  are 
weakest  over  the  Kuroshio  and  strongest  in  the  Southern 
Hemisphere  regions,  as  quantified  by  the  cross-correlation 
coefficients  between  the  crosswind  and  downwind  di¬ 
rection  gradients  and  the  crosswind  and  downwind  SST 
gradients,  respectively,  in  Table  3.  Additionally,  the 
magnitudes  of  the  correlation  coefficients  between  the 
downwind  direction  gradients  and  the  downwind  SST 
gradients  are  consistently  smaller  than  those  between 
the  crosswind  direction  gradients  and  the  crosswind  SST 
gradients.  However,  both  are  much  smaller  than  the 
correlations  between  the  crosswind  and  downwind  gra¬ 
dients  of  wind  speed  and  SST,  thus  explaining  the  higher 
variability  of  the  vorticity  and  divergence  binned  aver¬ 
ages  in  Fig.  11  compared  with  the  speed-only  contribu¬ 
tions  to  the  vorticity  and  divergence  in  Fig.  8. 

The  dependencies  of  the  direction  gradient  terms  on 
O'  quantified  statistically  in  Fig.  14  confirm  that  the 
crosswind  and  downwind  direction  gradients  vary  as  the 
sine  and  cosine  of  O',  respectively.  As  expected  in  all 
four  regions,  the  downwind  direction  gradient  response 
to  O'  (gray  curve  in  Fig.  14)  shows  that  the  surface  wind 
tends  to  rotate  anticyclonically  when  the  large-scale  flow 
is  from  cool  to  warm  water  (i.e.,  when  O'  =  0°)  and  cy¬ 
clonically  when  the  large-scale  flow  is  from  warm  to  cool 
water  (i.e.,  when  O'  =  ±180°).  The  crosswind  direction 
gradient  response  to  O'  indicates  a  maximum  surface 
diffluent  tendency  as  winds  blow  approximately  along 
perturbation  SST  isotherms  (i.e.,  when  O'  =  ±90°). 

Statistically,  the  direction  gradient  terms  can  thus  be 
represented  by 

=afM^cos(e'-<Af)  and  (11) 

(12) 

where  and  are  coupling  coefficients  for  the 
downwind  and  crosswind  direction  gradients,  respectively, 
and  and  are  phase  shifts  (see  Table  4).  These 
statistical  relations,  along  with  the  statistical  relations  for 
the  speed  gradients  in  Eqs.  (9)  and  (10),  are  used  in  sec¬ 
tion  4c  to  aid  in  understanding  the  differences  between 
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Fig.  10.  Binned  scatterplots  of  the  perturbation  crosswind  (black)  and  downwind  (gray)  speed  gradients  as 
functions  otO' .  The  bin  averages  were  computed  from  monthly- averaged  QuikSCAT  wind  and  AMSR-E  SST  fields 
over  the  75-month  period  considered  here  using  all  available  data  points  in  the  regions  enclosed  in  Figs.  1  and  2.  The 
points  and  error  bars  represent  the  means  and  standard  deviations  within  each  bin,  respectively.  The  solid  curves  are 
least  squares  fits  of  the  points  to  phase-shifted  sinusoids. 


the  vorticity  and  divergence  responses  to  the  crosswind 
and  downwind  SST  gradients. 

c.  Why  is  the  vorticity  response  to  crosswind  SST 
gradients  weaker  than  the  divergence  response  to 
downwind  SST  gradients? 

The  statistical  relations  developed  above  provide  an 
explanation  of  why  the  vorticity  response  to  crosswind 
SST  gradients  is  weaker  than  the  divergence  response  to 
downwind  SST  gradients.  Recall  that  the  vorticity  and 
divergence  dependencies  on  the  crosswind  and  down¬ 
wind  SST  gradients  were  interpreted  conceptually  from 


consideration  of  the  effects  of  SST  on  wind  speed  only 
(Chelton  et  al.  2001,  and  our  subsequent  studies).  From  a 
qualitative  perspective,  the  crosswind  direction  and  down¬ 
wind  speed  gradient  responses  to  6'  (the  black  and  gray 
curves  in  Figs.  14  and  10,  respectively)  are  nearly  in  phase, 
yielding  a  larger  amplitude  divergence  response  to  6'  than 
that  of  the  downwind  speed  gradient  alone;  in  contrast,  the 
downwind  direction  and  negative  crosswind  speed  gradient 
responses  to  6 '  (the  gray  and  black  curves  in  Figs.  14  and  10, 
respectively)  are  nearly  in  antiphase,  yielding  a  smaller- 
amplitude  vorticity  response  to  6'  than  that  of  the  cross- 
wind  speed  gradient  alone.  The  SST-induced  direction 


272 


Volume  23 


JOURNAL  OF  CLIMATE 


Table  4.  Values  of  the  coupling  coefficients  and  phase  angles 
appearing  in  Eqs.  (13)-(18)  computed  from  the  QuikSCAT  wind 
and  AMSR-E  SST  fields.  The  units  of  the  coupling  coefficients  are 
in  m  s~^ 
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gradient  perturbations  thus  enhance  the  divergence  re¬ 
sponse  to  downwind  SST  gradients  while  simultaneously 
reducing  the  vorticity  response  to  crosswind  SST  gradients. 

This  result  can  be  shown  more  rigorously  by  noting 
that  the  crosswind  and  downwind  speed  gradient  re¬ 
sponses  to  O'  (Fig.  10)  expressed  in  Eqs.  (9)  and  (10)  and 
the  crosswind  and  downwind  direction  gradient  responses 
to  O'  (Fig.  14)  expressed  in  Eqs.  (11)  and  (12)  can  be 
combined  to  relate  the  vorticity  and  divergence  fields  to 
Mt  and  O'  as 

V  X  u'  =  sm(e'  + 

V 

Crosswind  speed  gradient 

+  afM^cos{e'  -  4>f)  , 

Downwind  direction  gradient 

V  •  u'  =  af^Mj  cos(0'  + 

Downwind  speed  gradient 

+  af^Mj-sinie'  +<pf). 

Crosswind  direction  gradient 

The  components  of  the  vorticity  and  divergence  that 
each  term  represents  empirically  in  these  equations  are 
indicated  by  the  underbraces.  Using  trigonometric  iden¬ 
tities  and  rewriting  these  equations  in  terms  of  the  cross- 
wind  and  downwind  SST  gradients  using  Eq.  (7)  yields 

vx„'.„r(0+«r(0  and 


and  (13) 


(14) 


V  • 


where  the  set  of  coupling  coefficients  and 

are  defined  as 


^vort  ^  _^spd  ^Qg^spd  +  ^dir 

„vort  ^  _^spd  +  ^dir  (16) 

^  and  (17) 

af^  =  -aj^  sine/) +  af  cos(/)f .  (18) 

Consider  Eq.  (15)  for  which  is  the  coupling  co¬ 
efficient  for  the  vorticity  response  to  the  crosswind  SST 
gradient.  Both  and  sin(/)J^  are  statistically  found  to 
be  positive  quantities  (Table  4)  over  all  four  regions. 
Therefore,  the  downwind  direction  gradients  reduce  the 
vorticity  response  to  crosswind  SST  gradients  by  a  factor 
of  sin(/)J^  Likewise,  from  consideration  of  Eq.  (17) 
for  the  crosswind  direction  gradients  enhance  the 
divergence  response  to  downwind  SST  gradients  by  a 
factor  of  sin<pf^  since  and  sin<pf^  are  both  found 
statistically  to  be  positive  quantities  (Table  4)  over  all 
four  regions.  SST-induced  wind  direction  gradients  thus 
simultaneously  reduce  the  vorticity  response  to  cross- 
wind  SST  gradients  and  enhance  the  divergence  re¬ 
sponse  to  downwind  SST  gradients. 

It  is  also  noted  that  in  coupled  ocean-atmosphere 
simulations  of  tropical  instability  waves,  Seo  et  al.  (2007) 
and  Small  et  al.  (2009)  have  found  that  Ekman  pumping 
anomalies  associated  with  SST-induced  wind  stress  curl 
perturbations  tend  to  reduce  the  wind  stress  curl  relative 
to  the  wind  stress  divergence  perturbations.  This  ap¬ 
parently  occurs  because  surface  ocean  currents  tend  to 
be  quasi-nondivergent,  thus  affecting  the  wind  stress 
curl  more  than  the  wind  stress  divergence.  It  is  currently 
unknown  whether  similar  ocean  feedbacks  act  over 
midlatitudes  to  reduce  the  vorticity  response  to  cross- 
wind  SST  gradients  in  addition  to  the  wind  direction 
effects  reported  here. 


6.  Conclusions 

Analysis  of  the  differences  between  the  surface  vortic¬ 
ity  and  divergence  responses  to  SST  have  revealed  a  sig¬ 
nificant  response  of  the  wind  direction  to  mesoscale  SST 
perturbations  in  addition  to  the  well-known  positive  cor¬ 
relation  and  linear  response  of  the  surface  wind  speed  to 
SST.  Analysis  of  75  months  of  QuikSCAT  surface  vector 
winds  and  AMSR-E  SST  observations  showed  that  surface 
winds  accelerate  and  turn  anticyclonically  when  the  surface 
flow  is  from  cool  to  warm  SST  and  decelerate  and  turn 
cyclonically  when  the  surface  flow  is  from  warm  to  cool 
SST,  with  characteristic  wind  speed  changes  of  1-2  m  s~^ 
and  wind  direction  changes  of  4°-8°.  SST-induced  wind 
direction  perturbations  are  not  collocated  with  SST  per¬ 
turbations,  as  is  the  case  with  the  wind  speed  perturbations. 
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Fig.  11.  As  in  Fig.  3,  but  for  the  QuikSCAT  perturbation  vorticity  as  a  function  of  the  negative  perturbation 
crosswind  SST  gradient  -{dTIdn)'  (black)  and  the  QuikSCAT  perturbation  divergence  as  a  function  of  the  per¬ 
turbation  downwind  SST  gradient  (dTIds)'  (gray). 


but  are  displaced  poleward  by  about  1°  latitude.  The 
influence  of  mesoscale  SST  perturbations  on  surface 
wind  direction  has  not  been  previously  noted. 

The  vorticity  and  divergence  perturbations  are  found 
to  be  related  linearly  to  the  crosswind  and  downwind 
components  of  the  SST  gradient,  respectively,  as  in 
previous  studies.  The  vorticity  response  to  crosswind 
SST  gradients  is  30%-50%  weaker  than  the  divergence 
response  to  downwind  SST  gradients.  Although  this  is 
the  first  study  to  consider  the  responses  of  the  vorticity 
and  divergence  to  SST,  the  results  here  are  consistent 
with  previous  observational  studies  based  on  wind  stress 
curl  and  divergence,  as  discussed  in  the  introduction.  By 
explicitly  relating  the  vorticity  and  divergence  to  cross¬ 


wind  and  downwind  gradients  of  wind  speed  and  direction 
using  natural  coordinates,  we  showed  that  although  wind 
speed  gradients  contribute  nearly  equally  to  the  vorticity 
and  divergence  responses  to  SST,  the  SST-induced  wind 
direction  contributions  are  diametrically  opposed;  SST- 
induced  crosswind  and  downwind  gradients  in  wind  di¬ 
rection  reduce  the  vorticity  response  to  crosswind  SST 
gradients  while  simultaneously  enhancing  the  divergence 
response  to  downwind  SST  gradients. 

The  observational  results  shown  here  go  a  long  way 
toward  unifying  the  various  statistical  relationships  ob¬ 
served  among  mesoscale  perturbations  of  SST,  wind 
speed,  wind  direction,  vorticity,  and  divergence  fields.  At 
the  heart  of  these  interrelationships  are  the  SST-induced 
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Fig.  12.  As  in  Fig.  10,  but  for  the  QuikSCAT  perturbation  vorticity  (black)  and  divergence  (gray).  The  phase 
shifts  relative  to  pure  sine  and  cosine  curves,  respectively,  are  listed  in  the  lower  right  of  each  panel,  with  units  of 
degrees. 


responses  of  the  wind  speed  and  direction  to  mesoscale 
SST  perturbations.  The  relationships  expressed  here  in 
terms  of  the  surface  wind  speed  and  derivative  wind 
fields  have  been  extended  to  SST-induced  perturbations 
of  the  corresponding  wind  stress  fields  (OCE).  Because 
of  the  nonlinearity  of  the  stress-wind  relationship,  spa¬ 
tially  filtered  perturbations  of  the  wind  stress  curl,  wind 
stress  divergence,  and  wind  stress  magnitude  are  related 
to  those  of  the  vorticity,  divergence,  and  wind  speed,  re¬ 
spectively,  but  are  modulated  by  the  strength  of  the  am¬ 
bient  large-scale  wind  speed. 

The  feedback  onto  the  ocean  of  the  SST-induced  wind 
stress  curl  and  upper  ocean  mixing  is  the  focus  of  ongoing 
research.  It  is  becoming  apparent  that  ocean  circulation  is 


significantly  impacted  by  the  presence  of  SST-induced 
mesoscale  wind  variability  (e.g.,  Milliff  et  al.  1996;  Seo 
et  al.  2007;  Spall  2007a;  Jin  et  al.  2009;  Hogg  et  al.  2009; 
Small  et  al.  2009).  The  statistical  relations  among  vortic¬ 
ity,  divergence,  and  SST  developed  here  are  being  used 
to  develop  statistical  parameterizations  of  the  mesoscale 
surface  wind  stress  response  to  temporally  evolving  SST 
fields  in  ocean  models  to  simulate  the  effects  of  two-way 
coupling  between  the  ocean  and  atmosphere.  This  method 
may  be  used  to  simulate  mesoscale  wind  stress  vari¬ 
ability  in  long  time  series  of  surface  wind  stress  analyses 
from  numerical  weather  prediction  models,  which  un¬ 
derestimate  the  SST  influence  on  wind  stress  by  about 
a  factor  of  2  (Maloney  and  Chelton  2006;  Song  et  al.  2009). 
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Fig.  13.  (top  maps  in  each  panel)  The  perturbation  crosswind  direction  gradient  (colors)  overlaid  with  contours  of 
the  crosswind  SST  gradient,  (bottom  maps  in  each  panel)  The  perturbation  downwind  direction  gradient  (colors) 
overlaid  with  contours  of  the  downwind  SST  gradient.  The  data  in  each  map  were  averaged  over  the  75-month  period 
June  2002-August  2008.  The  contour  interval  for  the  perturbation  crosswind  and  downwind  SST  gradients  is  0.2°C 
per  100  km,  and  the  zero  contour  has  been  omitted  for  clarity.  Solid  and  dashed  contours  correspond  to  positive  and 
negative  values  of  the  SST  gradient,  respectively. 
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APPENDIX 

Physical  Interpretation  of  0' 

The  angle  6'  defined  by  Eq.  (8)  is  important  in  quan¬ 
tifying  the  response  of  the  surface  perturbation  vorticity 
and  divergence  fields  to  spatially  varying  SST.  Non- 
linearities  and  spatial  high-pass  filtering  of  the  crosswind 
and  downwind  SST  gradient  fields  make  exact  analytical 
formulation  of  6'  in  terms  of  wind  and  SST  intractable. 
Nevertheless,  an  effort  is  made  here  to  clarify  the  prac¬ 
tical  physical  relevance  of  6'  to  the  orientation  of  the 
surface  winds  to  SST  fronts,  which  conceptually  control 
spatial  gradients  in  wind  speed  and  direction  as  discussed 
throughout  this  study. 
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Fig.  14.  As  in  Fig.  10,  but  for  the  perturbation  crosswind  (black)  and  downwind  (gray)  direction  gradients  as 

functions  of  6'. 


The  unfiltered  crosswind  and  downwind  SST  gradient 
components  can  be  expanded  as  the  sum  of  spatially 
high-pass  filtered  components,  represented  by  primes, 
and  smoothed,  spatially  low-pass  filtered  components, 
represented  by  tildes: 


dn  \dn  J 


and 


(Al) 


(gj.ivrisinr, 


(A3) 

(A4) 

(A5) 


dT 

l)s 


(A2) 


|Vr|sin6l', 


(A6) 


An  empirical  approximation  can  then  be  made  by  writing 
the  crosswind  and  downwind  SST  gradient  components 
in  the  following  forms; 


IVri  sind,  and 


(A7) 


IVT'lsin(e')  (°C  per  1 00km)  IVT'lcos(e')  {°C  per  1 00km) 


15  January  2010 


O’NEILL  ET  AL. 


277 


3 

2 

1 

0 

-1 

-2 

-3 

-3-2-10123 
Perturbation  Downwind  SST  Gradient 
(°C  per  100km) 


a) 


1  1  1  1  1 

/ 

/■ 

1_ 1_ L 

1 

b) 


Smoothed  Downwind  SST  Gradient 


(°C  per  100km) 


c) 


Perturbation  Crosswind  SST  Gradient 


(°C  per  100km) 


d) 


Smoothed  Crosswind  SST  Gradient 
(°C  per  100km) 


Fig.  A1.  Binned  scatterplots  of  the  relations  expressed  by  Eqs.  (A5)-(A8).  (a)  |Vr'|  cosO'  vs.  (dTIds)';  (b) 
|Vr|cos0vs.  (dT/ds);  (c)  |Vr'|  sin0  vs.  (aT/a^)';  (d)  |Vr|  sina  vs.  (aT/a^z).  The  points  and  error  bars  in  each  panel 
represent  the  means  and  ±  1  std  dev  within  each  bin  computed  from  the  first  monthly  averaged  AMSR-E  SST  and 
QuikSCAT  wind  fields  used  in  this  study  over  the  Agulhas  Return  Current.  The  dashed  line  in  each  panel  has  unit 
slope  and  is  shown  for  reference. 


®  =  |V7’|cos0, 

(A8) 

9  =  tan  ^ 

(dTIdn) 

ydsj 

(df/fs) 

where  6  is  defined  as  the  counterclockwise  angle  be-  Equations  (A5)-(A8)  are  only  approximations  that 
tween  the  surface  wind  vector  and  \T,6'  is  defined  as  in  cannot  be  rigorously  justified  analytically  because  the 
Eq.  (8),  and  6  is  defined  analogously  as  loess  smoother  used  to  produce  the  spatially  filtered  SST 
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VT' 


Fig.  A2.  Schematic  of  the  orientation  of  the  surface  wind  vector  u,  the  SST  gradient  vector 
Vr,  and  the  perturbation  SST  gradient  vector  VT'.  The  counterclockwise  angle  6'  -  0  is  as 
deduced  from  Eq.  (A16). 


fields  here  is  not  commutative  with  the  nonlinear  vector 
magnitude  operator.  However,  they  are  shown  empiri¬ 
cally  in  Fig.  A1  to  be  nearly  exact,  despite  the  non- 
linearities,  based  on  terms  computed  from  the  AMSR-E 
SST  and  QuikSCAT  wind  fields  averaged  over  the  1-week 
period  29  December  2002-4  January  2003,  which  is  rep¬ 
resentative  of  the  6+ -year  data  record  used  here. 

An  interpretation  of  6'  can  be  obtained  by  eliminating 
9  and  T  from  consideration  and  solving  for  6’  in  terms  of 
0  in  a  form  that  reveals  that  6'  is  the  angle  between  the 
wind  vector  u  and  the  perturbation  SST  gradient  vector 
\T.  With  the  above  forms,  Eqs.  (Al)  and  (A2)  can  be 
written  as 

|Vr|  sin^  =  |Vr'|  sin^'  +  |Vr|  sinO  and  (A9) 

I  vr|  cose  =  \\r  \  cose'  +  \\f\  cose.  (aio) 
Squaring  and  adding  these  two  equations  yields 
\\T'f  +  2|V7’'||vf|  cos(e'  -e)+  |vU  =  ivrp- 
Hence, 

(All) 


where 

|vrf-|vrp-|vff 

2|vr||vr|  ‘  ^ 

Substituting  Eq.  (All)  into  Eqs.  (A9)  and  (AIO)  with 
the  trigonometric  identity  sin(cos“^R)  =  ±(1  - 
and  summing  sine  and  cosine  trigonometric  identities 
yields 

|Vr|sin0=|Vr|sin0' 

+  |Vr|[R  sin^'  ±  (1  cos^']  and 

\\T\ COS0  =  I  VT'  I  cosd'  +  \\f\[Rcose'  ±  (1  -  sine']. 

These  can  be  rewritten  as 

\\ T\  sine  =  A  sine' ±  B  cose'  and  (A13) 
|Vr|  cos^  =  Acos^' ±  5sin^',  (A14) 

where  A  =  (|Vr'|  +  \\f\R)  and  B  =  \\f\(l-R^f^. 
The  positive  or  negative  root  in  each  equation  means 
there  are  four  possible  sets  of  equations.  However,  by 


9  =  0'-  cos  ^R, 
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Fig.  A3,  (top)  Binned  scatterplot  of  the  angle  6'  between  the 
surface  wind  vector  and  the  perturbation  SST  gradient  vector 
binned  as  a  function  of  0'.  (bottom)  Histogram  of  the  differences 
between  O'  and  O';  the  distribution  has  a  standard  deviation  of  7.6°. 
The  points  and  error  bars  in  the  top  panel  represent  the  means  and 
±1  standard  deviation  within  each  bin  computed  from  the  first 
monthly-averaged  AMSR-E  SST  and  QuikSCAT  wind  fields  used 
in  this  study  over  the  Agulhas  Return  Current.  The  dashed  line  in 
the  top  panel  has  unit  slope  and  is  shown  for  reference. 


noting  that  =  |Vr|,  only  two  sets  of  equations 

are  possible,  namely 

\\T\sm9  =  AsinO' B  cos9'  and 

|Vr|  cos9  =  A  cos9'  -  B  sm9', 


)  '  \  dy  J  ' 


\^Tr=  —  +  — 


ary 


and 


I'dT 

\ax  dx  J  \dy 


Using  these  relations,  and  after  some  manipulation 
and  simplification,  Eq.  (A15)  can  be  written  in  the  com¬ 
pact  form 


\T  •  vr' 

cos(6i  -e)- 


(A16) 


From  the  definition  of  the  dot  product  in  the  numer¬ 
ator  of  Eq.  (A16),  it  is  evident  that  the  angle  9'  —  9  is  the 
counterclockwise  angle  between  the  vectors  \T  and 
Vr.  A  schematic  of  these  angles  and  SST  gradient 
vectors  is  shown  in  Fig.  A2.  From  this,  it  can  be  deduced 
that  9'  is  the  counterclockwise  angle  between  u  and  the 
perturbation  SST  gradient  vector  \T.  This  deduction  is 
confirmed  statistically  in  the  top  panel  of  Fig.  A3,  where 
the  angle  between  u  and  denoted  as  is  binned  as 
a  function  of  the  angle  9'  computed  from  Eq.  (8)  using 
the  6+  years  of  monthly-averaged  QuikSCAT  wind  and 
AMSR-E  SST  fields  utilized  in  this  study.  Differences 
are  quantified  in  the  histogram  in  the  bottom  panel  of 
Fig.  A3,  which  shows  a  very  narrow  distribution  cen¬ 
tered  at  0°  with  a  standard  deviation  of  7.6°.  It  is  thus 
concluded  that  the  angle  between  the  surface  wind  and 
perturbation  SST  gradient  vectors  is  very  well  repre¬ 
sented  by  the  angle  9'. 


or 
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